APPROXIMATION BY POLYNOMIALS AND BLASCHKE PRODUCTS 
HAVING ALL ZEROS ON A CIRCLE 



DAVID W. FARMER AND PAMELA GORKIN 

Abstract. We show that a nonvanishing analytic function on a domain in the unit disc 
can be approximated by (a scalar multiple of) a Blaschke product whose zeros lie on a 
prescribed circle enclosing the domain. We also give a new proof of the analogous classical 
result for polynomials. A connection is made to universality results for the Riemann zeta 
function. 



1. Introduction 

While every analytic function on a domain can be approximated pointwise on the domain 
by a polynomial, there are many other interesting questions that can be asked about the 
approximating polynomials; for example, can the sequence of approximating polynomials be 
chosen to be uniformly bounded? In this paper, we focus on the case in which the domain 
is the open unit disc and we ask where the zeros of the polynomials will lie. 

When our attention is focused on bounded analytic functions, there is another class of 
functions that can be used to approximate our bounded analytic function: the set of Blaschke 
products. A finite Blaschke product is a function of the form 

N 

3=1 J 

where \a,j\ < 1 for all j and |A| = 1. Caratheodory's theorem (see [7J p. 6], for example) 
shows that if / is an analytic function defined on the open unit disc, D, and / is bounded 
by 1 in modulus, then there is a sequence {Bk} of finite Blaschke products converging to / 
pointwise on D. Again, it is certainly interesting to ask where the zeros of the approximating 
Blaschke products may lie. 

Looking at more general domains, one natural question is the following: Given a holomor- 
phic function in a Jordan region, when can it be approximated by a polynomial with zeros 
lying on the boundary? This question was answered by G. MacLane [16J in 1949. Curiously, 
many texts dealing with the study of polynomials or approximation by polynomials do not 
include reference to this work ([17], [19]). though three different proofs of this result are 
given by Korevaar [T2] . 
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MacLane's work focused on showing that a zero-free holomorphic function can be approx- 
imated by a polynomial with zeros on the boundary, when the boundary satisfies certain 
smoothness conditions. C. Chui [2], [3] looked at the problem of bounded approximation of 
a zero-free bounded holomorphic function by what he called C-polynomials. Chui showed 
that every zero-free bounded holomorphic function defined on D, can be boundedly approx- 
imated by polynomials with zeros lying on the unit circle. In 1968, Z. Rubinstein showed 
that given a zero- free holomorphic function with /(0) = 1, there exists a sequence of C- 
polynomials mapping to 1 that converges to / uniformly on every compact set in D. In 
addition, when the function / is bounded, the sequence converges boundedly. One natural 
approach, that of looking at the zeros of the partial sums of a series, has been further studied 
by Korevaar and others. This and related work can be found in [T3], P3], and |14j . 

In this paper, we will focus on polynomials as well as Blaschke products. In general, a 
Blaschke product is a function of the form 

oo 

where n is a nonnegative integer, < |eij| < 1, |A| = 1, and YljO- ~ \ a j\) < 00 • As indicated 
above, when there are finitely many aj, the function is said to be a finite Blaschke product; 
the term Blaschke product is used when the distinction between finite and infinite products 
is unimportant. The Blaschke product is determined by its zeros and every bounded analytic 
function on D can be written as a product of a Blaschke product and an analytic function 
with no zeros on D. In fact, every bounded analytic function / has what is known as an 
inner-outer factorization; that is, there exists an analytic function J, bounded by 1 in the 
unit disc, with |/| = 1 almost everywhere on the unit circle, and an outer function g with the 
property that / = Ig. The inner function I can be further factored into a (possibly infinite) 
Blaschke product, defined with the zeros of the function /, and a singular inner function. 
Though the singular inner function will have no zeros in D, Frostman's theorem [7J p. 79], 
shows that it can be uniformly approximated by a Blaschke product. Therefore, every inner 
function can be uniformly approximated by Blaschke products. 

R. Douglas and W. Rudin [2] showed that such results can be applied to another very 
important algebra of functions. In fact, as they showed, every essentially bounded measurable 
function u with modulus 1 a.e. on the unit circle can be approximated using Blaschke 
products and their conjugates. More precisely, given e > 0, there exist Blaschke products 
b\ and 62 such that \\u — &1&2II < £ - Several interesting results on these algebras followed, 
and then, in 1982, P. Jones [10] showed that the Blaschke products above could be chosen to 
be interpolating Blaschke products; that is, the zero sequence {z n } of the Blaschke product 
has the property that given any bounded sequence {w n } of complex numbers, there exists 
a bounded analytic function / such that f(z n ) = w n . A characterization of such sequences 
due to L. Carleson pQ makes such sequences quite easy to deal with. In particular, the 
zeros of the interpolating Blaschke products are separated in an extremely useful way, and 
as a consequence interpolating Blaschke products have many desirable properties. In trying 
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to understand these results, Jones and Garnett [3, p. 430] asked whether every Blaschke 
product can be uniformly approximated by interpolating Blaschke products. While many 
interested related results exist (see, for example, [8] and [9] for recent results ) this problem 
remains open. Our interest in controlling the zeros of approximating Blaschke products is 
directly related to this very important question. If we can control the placement of the zeros 
of an approximating Blaschke product, we have a chance of constructing approximating 
interpolating Blaschke products. 

Our main result (Corollary 12. 3p is that given an analytic function g that has no zeros in a 
neighborhood of {z : \z\ < r}, for all e > and 5 > there exists a constant cb and a finite 
Blaschke product B with all zeros on the circle {z : \z\ = r} such that 

\g(z) - c B B(z)\ <e (1.3) 

on {z : \z\ < r — 5}. Our approach also provides a new and relatively simple proof of the 
fact that a nonvanishing analytic function can be approximated uniformly on compacta by 
polynomials having zeros on a prescribed circle. 

The paper is organized as follows. In Section [2] we consider the special case of discs 
centered at 0. Then in Section l37fl we use the special case to prove more general cases. In 
Section H] we discuss the relationship with universality results for the Riemann zeta function 
and its connection to random matrix theory. 



2. Approximation around 

Here is a very simple proof that a polynomial p that does not vanish on a neighborhood 
of the closure of the unit disc can be approximated uniformly on compacta by polynomials 
with zeros on the unit circle: 

Note that if the degree of p is m, then the polynomial p*(z) = z m p(l/z) has all zeros 
inside the unit disc. Let 

B(z)=p*(z)/p{z). (2.1) 

Then B is analytic in D, continuous on the unit circle, maps D to itself, the unit circle to 
itself, and the complement of the closed disc to itself. Therefore B is a Blaschke product. 
Now for fceN, the set of points in D for which B{z) = z~ k lie on the unit circle. Therefore, 
the polynomial 

p(z) + z k p*(z) (2.2) 

has all its zeros on the unit circle and approximates p on compacta as k — > oo. This is, more 
or less, the proof given by Z. Rubinstein [18j in 1968. However, it does not seem possible 
to adapt this proof to the case of Blaschke products. We present an alternate proof of this 
result, before turning to our theorem on Blaschke products. 

First we consider the case of approximating on a disc centered at 0. 
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Theorem 2.1. Suppose f is analytic and nonvanishing in a neighborhood of \z\ < r. Then 
there exist numbers = \r]j\ — 1, A G C ; and positive integers such that 

oo 

/( 2 )=A[J(l+(/)^(l + , ) /), (2.3) 
i=i 

for \z\ < min{r, 1}, with the convergence uniform on \z\ < min{r, 1 — 5} for any 5 > 0. 

In particular, f can be approximated on \z\ < min{r, 1 — 5} by polynomials having all roots 
on the unit circle. 

Theorem 2.2. Under the same conditions as Theorem \2.1\ if R < 1 then 

v{j) 1 + Vj zi 



(2.4) 



for \z\ < min{r, 1} ; with the convergence uniform on \z\ < min{r, 1 — 5} for any 5 > 0. 

Corollary 2.3. Suppose g is analytic and nonvanishing in a neighborhood of \z\ < r < 1. 
For all e > and 5 > there exists a constant cb and a Blaschke product B having all zeros 
on \z\ = r such that 

\g(z) - c B B(z)\ < e (2.5) 

for \z\ < r — 5. 

Before giving the proof of Theorems 12.11 and 12.21 we describe the construction in a context 
that avoids the issue of convergence. 

2.1. Formal power series as infinite products. 

Proposition 2.4. Suppose 

oo 

f( z ) = j2 a - zn ( 2 - 6 ) 

n=0 

is a formal power series with a n e C and a ^ 0. We can write 

oo 

/(*) = «o I^ 1 + ^r 0) (l + (2-7) 

where |^| = j^l = 1 and z/(j) is a nonnegative integer. 
Proposition 2.5. Suppose 

oo 

/(z) = £a n z n (2.8) 

n=0 

a formal power series with a„6C and a 7^ 0, an<i let R < 1. We can write 

/M-ofif 1+ ^.r 1 + . (^) 
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where = \r]j\ = 1 and u(J) is a non-negative integer. 

Corollary 2.6. Suppose g is analytic in a neighborhood of with g(0) ^ 0, and suppose 
< r < 1. Then for all J > there exists a constant cj and a Blaschke product Bj having 
all zeros on \z\ = r such that 

g(z)-cjBj(z)=0(z J ), (2.10) 

as z — > 0. 

The proof of Corollary 12.61 is the same as the proof of Corollary 12.31 The proof of latter 
will be presented in Section I2.3[ following the proof of Theorem 12.21 

We define the product representations inductively, using the following Lemma. 

Lemma 2.7. Let R Q > 0. Every w G C can be written as w = m£ + rj where |£| = 
\tj\ = Rq and m is a positive integer. There are four such representations if \w\ > Rq, two 
representations if < \w\ < Rq, and infinitely many ifw = 0. 

Proof. We need only consider the case R = 1. 

Consider the unit circle centered at m e 10 . As 9 goes from to 2tt, that circle sweeps out an 
annulus whose inner and outer radii differ by 2. Every point in the interior of that annulus 
has two representations of the form m e ld + 77. As m varies, those annuli cover the complex 
plane. Every point lies in the interior of two of the annuli, except the points \w\ < Rq or 
whose absolute value is an integer. □ 

Proof of Proposition \2.J\ We may assume = 1. By Lemma 12.71 with Rq = 1, we can 
choose £1, T)i, and u(l) so that a\ = + rji. 

Set 

P l {z) = {l+i 1 zY^{l + r ll z) (2.11) 
and note that P\{z) = 1 + a^z + 0(z 2 ). Therefore, 

(2 - 12) 

where b = 1 and bi = 0. Now choose £ 2 , V2, and z/(2) so that b 2 = v{2)^ 2 + V2- Setting 

P 2 (z) = (l + t 2 z 2 yM(l + r l2 z 2 ) 

= l + b 2 z 2 + 0(z 3 ) (2.13) 

we have 

Pi{z)P 2 {z) = 1 + a lZ + a 2 z 2 + 0{z 3 ). (2.14) 
Proceeding inductively we obtain f(z) = Y\jPj(z). □ 

Proof of Proposition \2.5[ Note that 



1 + ijWzi 



1 + (1-R j )^z j + 0(z 2j ). (2.15) 
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So everything goes through in the previous proof, with the modification that we apply 
Lemma O with R Q = I - W □ 

The proofs given above are just formal calculations, and it is not clear what convergence 
properties the infinite products might have. Even if / represents an analytic function in a 
neighborhood of the origin, the products can only converge where / does not vanish. The 
convergence of the products will depend on the growth of the numbers and the above 
constructions do not appear to shed light on this. 

In the next section we organize the proof in a different way that may appear more cum- 
bersome, but it gives information about the analytic properties of the infinite product. 

2.2. Proof of Theorem HZQ 

Proof of Theorem \2.1\ Let 

J -(z) = J2 a nZ n - (2-16) 

J n=0 

Since /'// is analytic in a disc slightly larger than \z\ < r, there is a C > and k > 1/r so 
that 

\a n \ < Cn n (2.17) 

for all n. If r > 1 we set k = 1 + 5 for some 5 > 0. In particular, k > 1. 
Let 

j 

g{z) = gj(z) = + ^T 0) (l + Vi* 1 ), (2-18) 
i=i 

where = \rjj \ = 1 and u(j) is a non-negative integer. We will choose those parameters so 
that the first J terms in the Taylor series for g'/g match those of /'//. 
We have 



V (jiA.iX.y 1 1 — r + jmJ- 1 A 

oo > 

2- 1 ^ j ( ^(j)(-l) m+1 ef + V (m+1) + (-l) m+1 77f + V 
KKJ m=0 ^ 



fc=o i<j<J 

j|(fc+l) 



fc+i fc+i 



fe=0 



say. On the third line we changed summation index = j(m + 1) — 1, and the notation n|m, 
read "n divides m," means that m/n is an integer. 
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Now we show how to choose the parameters in g to match the Taylor series coefficients of 
the logarithmic derivatives. 

By Lemma 12.71 we can choose a non-negative integer u(l) and complex numbers = 
1 771 1 = 1 so that — a = f(l)£i + 771- Thus, b = a , and we have matched the first Taylor 
series coefficients of /'// and g'/g. 

Since the only positive integer that divides + 1 = 1 is 1, we will continue to have 60 — 
no matter what we later choose for rjj, and for j > 2. Likewise, once we have chosen 
£j, rjj, and u(J) for j < K so that aj = bj for j < K — 1, we will continue to have aj = bj for 
j < K — 1 because if j|(A; + 1) then j < k + 1. 

To choose £j, 77^, and ^(j) for j > 2, we have to deal with the fact that the j = 1 terms 
make a contribution to all of those coefficients. Similarly, the j = 2 terms contribute to all 
of the later even-index coefficients, and so on. Breaking the sum defining bx into two parts, 
we find 



>K 



K+l / K+1 K+1 \ K+1 / K+1 K+1 \ 



i<j<K x x k+i<?'<J 

i |(X+i) j l (K+i) 



K+1 K+1 



= - 2. ( n ! ./ K/K, ' +v ) + + 1) + + w+ij • 

i|(Jf+i) 

(2.20) 

The terms in the sum on the second line of (12.201) have already been chosen, so we can use 
Lemma to choose v{K + 1), £#+1, and r/x+i so that bx = clk- 

Proceeding in this way we match the first J coefficients of the logarithmic derivatives. 

It remains to bound v{j) so that we can bound the tail of ( 12. 19ft . 

By (I2.20p and the fact that bx = clk we have 

(K + l)u(K + l)<K + l + \a K \+ J2 (2-21) 

•7'KK+l) 
l<i<K 



By Lemma 12.81 given at the end of this section, the above estimate implies that there 
exists C so that nvin) < C ' k u for all n > 1. This is sufficient to estimate the tail for 
\z\ < 1/k because the coefficient of z n in ( I2.19P is bounded by 

IU<^l>(j)+j)<< jV > (2-22) 

3<J 

where we use <C to mean the quantity is bounded by a constant times J 2 k j . So 

|&n||z n |«JV \ z \ n < YZT-\ kJ \ z \ J , (2-23) 



n>J+l n>J+l 

which goes to as J — > 00 because \z\ < 1/k < 1. 
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This shows that g'(z)/g(z) is close to f'(z)/f(z) for \z\ < 1/k. We can antidifferentiate 
using Cauchy's theorem, so log(/) is close to log(gj) + cj for \z\ < 1/k,, for some constant cj. 
Now exponentiate to get that f(z) is close to e Cj g(z). Since g(0) = 1, choose cj = log(/(0)). 

This completes the proof of Theorem 12.11 □ 

Lemma 2.8. Suppose > for j > 0. If there exist k > 1 and C > so that 

(n + l)u(n + l) <n 2 + CK n + ^ (2.24) 

j\(n+l) 
\<j<n 

for all n > 1, then there exists C so that nv(n) < C'n n for all n > 1. 
Proof. First choose N so that 

" < - (2.25) 



K (n+l)/2 3 

if n > N. Then choose C so that 

(1) a > 3C, 

(2) n 2 < \C'K n for all n, and 

(3) n u(n) < C'K n for n < N. 

Note that ([2]) uses only the fact that k > 1, and ([3]) uses only that k > and N is finite. 

Now we prove the desired estimate by induction. Suppose n v{n) < C K n for n < M, 
where M > N, and suppose n — M + 1. Using the first two conditions on C, the induction 
hypothesis, and the fact that all proper divisors of n + 1 are at most (n + l)/2, we have 

{n + l)v{n + l)< l -C'K n + l -C'K n + MJ) 

j<(n+l)/2 

< -C'K n + V C'K j 

~ 3 ^ 

j<(n+l)/2 

~ 3 2 

2 cv* + cv +1 n + 1 



3 2k("+ 1 )/ 2 
< C'K n+1 . (2.26) 

The last inequality follows from n > N and the choice of N. That completes the proof of 
Lemma EE □ 



2.3. Proof of Theorem 12.21 We describe how to modify the the proof of Theorem 12.11 to 
give a proof of Theorem 12. 2[ 

Proof of Theorem \2.Si Let 
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h ' g 

where g is the function (12.181) appearing in the proof of Theorem 12.11 Writing 



we have 



ti 



c k 



k=0 



;i - R k+1 )b k , 



(2.28) 



(2.29) 



(2.30) 



where b k are the Taylor series coefficients of g'/g given in (I2.19p . 

Thus, when matching the coefficients of h'/h and /'//, everything goes as before if in each 
equation we replace b k by c k and a k by a k /{l — R k+1 ). So the choices of i] k , and u{k) 
follow the same steps. The final step of bounding v{K) involves replacing inequality (I2.2ip 
by 

(K + l)u(K + l) <K + 1 ! "' v 



1-R K - 



(2.31) 



j'KJf+i) 

1<3<K 



□ 



But that implies the bound we need on u(j) because the only fact we used about a k is 
\a k \ < Cn k for some C > 0. 
This completes the proof of Theorem 12.21 

We turn to the proof of Corollary 12.31 

Proof of Corollary \2.3[ In Theorem I2.2[ let R = r 2 and set g(z) = f(z/r), so 



1 + (jZ> 



sup 

\z\<i-s 



g(rz) = f(z) = a JJ ( - 

3=1 ^ ' Sr 

By Theorem 2.2, for each 5 > we may choose N so that 

g(rz) -a Yl ( — 



1 + 1]jZ j 

1 + 7]jr 2 i zi 

1 + T}jZ j 



+ ^r^izi J 1 + rjjr 2 i zi 



< e. 



It remains to rearrange the above product to recognize it as a Blaschke product. Letting 
w = rz we have 



sup 

|to|<r(l— 5) 



a 



n 



r)jr : 'w J 



< 6. 



Thus, 



sup 

\w\<r(l-6) 



N 



' ) n r 2j ^vj 



1 + tijrivji 



) 



1 + rjjr^w 3 



< e. 
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Letting otj = r J '£j and j3j = r^rjj, we have 



sup 

H<r(l-<5) 



N 



a 



< e, 



(2.32) 



where C(w) = TT ( ItL^Y® (i 



+ /W 



. Since each factor of C is a Mobius transfor- 



mation composed with u^, each factor is a Blaschke product and therefore C is a Blaschke 
product as well. □ 



3. Corollaries of Theorem 12.11 

We deduce some corollaries. 

3.1. Approximation on pseudohyperbolic discs. As a corollary to Theorem 12.11 we 
prove a result about approximating on other discs contained in the unit disc. We recall first 
that the pseudohyperbolic distance between two points z and w in D is defined to be the 
distance 

z — w 



p(z,w) 



1 



wz 



(3.1) 



For a G D and r with < r < 1 we let D p (a, r) = {z : p(a, z) < r}. Given a Euclidean disc 
D(cio,r ), we may rotate it so that the center lies on the positive real axis, and let x and y, 
with | x | < y, denote the points in which the bounding circle C intersects the real line. Let 
a a be the Mobius function 

' ' Z + a (3.2) 



OC„\Z) 



1 



az 



and let R = ±±g£. Then R > 1 and if a = R - VR 2 - 1 then r = -a'^x) = a~ l (y). 
Since a is real, a" 1 maps C onto a circle C\ passing through r and — r and since the real 
line is orthogonal to C, the real line must be orthogonal to C\. Therefore a^ 1 maps C onto 
{z : \z\ = r}. Thus, the disc D(a ,r ) is rotation of a pseudohyperbolic disc D p (a,r) for 
some a, r. This means that 



£>,(a,r) =a a (D(0,r)) 



(3.3) 



For basic information about automorphisms of the disc, see Garnett [7]. 

Theorem 3.1. Let f be a function that is analytic and nonvanishing in a neighborhood of 
the disc D(a ,r ) C D. Then f can be uniformly approximated on D(a ,r ) by a polynomial 
with zeros on the unit circle. 

Proof. Suppose / has no zeros in a neighborhood of the closure of a pseudohyperbolic disc 
D p (a,r). Then / o a a has no zeros in a neighborhood of the disc D(0,r). By Theorem 12. H 
there is a polynomial p with zeros on the unit circle such that 

\\foa a -p||d( ,p) < e. (3.4) 
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Therefore by (13. 3 j) and a change of variables, 

\\f -poa a ' l \\ DpM < e. (3.5) 
Now, letting z 1 , . . . ,z N denote the zeros of p on the unit circle we see that 

N 

poa-\z) =Y[(a a -\z) - Zj ). (3.6) 
i=i 

This is a rational function with poles outside the (closed) unit disc and zeros at a a (zj) for 
j = 1, . . . , N. Thus, the zeros of this rational function also lie on the unit circle. 

Now choose s < 1 so that D p (a, r) C D(0, s). Since poa' a 1 {z) is analytic and nonvanishing 
in a neighborhood of D(0, s), we can apply Theorem 12.11 again to get a polynomial q so that 

||g -po a a _1 || D ( , s ) < £, (3.7) 

which implies 

\\q-po a a ~ l \\ Dp (a,r) < e. (3.8) 
Combining (13. 5p and (13.81) gives ||/ — q\\D p {a,r) < 2e, as required. 

□ 



There is a Blaschke product version of this result that can be obtained in a similar, but 
simpler manner; that is, if we use Corollary 12.31 in place of Theorem 12.11 and note that Boa' 1 
is a finite Blaschke product (see [7J p. 6]) with zeros on the boundary of D p (a,r) whenever 
the zeros of B lie on {z : \z\ = r}, we obtain the following result. 

Theorem 3.2. Let f be a function that is analytic and nonvanishing in a neighborhood of 
the disc D(ao,r ) C D. Then f can be uniformly approximated on D(ao,r — d~) by a constant 
times a Blaschke product with zeros on the circle {z : \z — ao\ = To}- 

From these results, we obtain a corollary about functions with zeros. For < p < oo and 
/ an analytic function on D, we say that / £ H p if 

sup^ j \f(re^fde=\\f\\ p HP <oo. (3.9) 

It is well known that given a nonzero function / £ H p the zero sequence of /, denoted (z n ), 
is a Blaschke sequence. Letting C\ denote the (possibly infinite) Blaschke product with zeros 
(z n ) there exists a function g that is analytic on D and has no zeros in D such that / = C\g. 
Applying the previous theorem to g we obtain the following. 

Corollary 3.3. Let < p < oo and let f £ H p . If D(ao,ro) C D then f can be uniformly 
approximated on D(a , r —5) by functions of the form c CiC 2 where c is a constant, C\ is the 
Blaschke factor of f , and Ci is a Blaschke product with zeros on the circle {z : \z — a \ = r }. 
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4. The Riemann zeta function and random matrix theory 

We describe a new approach and another motivation for studying approximations by 
polynomials with zeros on the unit circle. This involves a combination of a universality 
result for the Riemann zeta function, and the connection between the zeta function and 
random matrices. 

We recall Voronin's universality result. 

Theorem 4.1. (Voronin [22J ED] ) Let < r < | and suppose g is a nonvanishing continuous 
function on the disc \s\ < r which is analytic in the interior. Then for any e > 0, 

liminf — meas It £ [0, T] : max |C(~ + %r + s) - g(s)\ < el > 0. (4.1) 

T^>oo T I \s\<r J 

Here ((s) is the Riemann zeta function, defined by 



oo 1 

coo = £^ (4-2) 



n 

n=l 

for dt(s) > 1. The standard reference is Titchmarsh EH and we merely cite a few facts 
that are relevant to our discussion here. The zeta function continues to a meromorphic 
function with a single simple pole at s — 1. The functional equation relates ((s) to £(1 — s). 
Combining the functional equation with (14. 2[) shows that the most interesting behavior of 
the zeta function will occur in the "critical strip" < 3?(s) < 1. Voronin's theorem says that 
the zeta function is universal in the right half of the critical strip. Furthermore, any given 
nonvanishing analytic function on a disc of radius < | is closely approximated by a positive 
proportion of shifts of £(| + is). 

The Riemann Hypothesis is the conjecture that the zeros of the zeta function in the 
critical strip all lie on the "critical line" dt(s) = \. Recently there has been significant 
progress in understanding (i.e., conjecturing) the behavior of the zeta function in the critical 
strip by modeling the zeta function by the characteristic polynomial of a random unitary 
matrix jTU Hj. Specifically, £(| + iT + z), for $t(z) > 0, is modeled statistically by A(e~ z ), 
where A is the characteristic polynomial of a random matrix from the unitary group U(N), 
chosen uniformly with respect to Haar measure, where N is approximately log(T/27r). 

Thus, Voronin's theorem and the above discussion suggests the theorem that characteristic 
polynomials of unitary matrices (ie, polynomials having all zeros on the unit circle) can 
approximate nonvanishing functions in the unit disc. But not quite: the normalization used 
in analytic number theory has A(0) = 1, corresponding to the fact that lirnsjj( s )_j. 00 ((s) = 1. 
Thus, characteristic polynomials can approximate nonvanishing functions on D(ao,ro) C D 
provided ^ D(a ,r ). This of course follows from the results described at the beginning of 
Section [2j but we suggest that it would be interesting to give a purely random matrix proof. 
That is, to understand the distribution of 

A' A( n ) \ 

_(*),...,_(*)! (4.3) 
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with < x < 1 fixed, for A the characteristic polynomial of a random matrix in U(N) 
as iV — > oo. The approximation theorems show that each z G C n is of the above form 
if N is sufficiently large. But a purely random matrix calculation, providing a probability 
distribution for (14. 3p . could give yet another proof. 

It would also be interesting to prove a random matrix analogue of the fact that a positive 
density of shifts of the zeta function approximate a given function. Suppose / is a nonvan- 
ishing analytic function on the unit disc with /(0) = 1, and suppose < r < 1 and e > 
are given. If e N matrices U G U(N) are chosen randomly with respect to Haar measure, the 
problem is to determine the probability that at least one of those e N matrices satisfies 

\det(I-U*z) -f(z)\ <e for all \z\ < r. (4.4) 

See j6] for an explanation of why e N matrices are chosen from N(U). Voronin's universality 
result suggests that if N is sufficiently large then the above probability is positive and 
bounded below independent of N. 

If this random matrix approach is successful, it will produce polynomials of a very different 
form than those constructed in the proof of Theorem 12.11 That proof involved polynomials 
with high multiplicity in their zeros. But the characteristic polynomial of a random unitary 
matrix has, with probability 1, only simple zeros, and those zeros tend to be very evenly 
spaced on the unit circle. 
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